Geometrical Properties of Loops and Cluster Boundaries 



John Cardy^ 

All Souls College 
and 

Department of Physics 
Theoretical Physics 
1 Keble Road 
Oxford 0X1 3NP, UK 

We discuss how the statistical properties of the area and radius of gyration of single 
self-avoiding loops, and of Ising and percolation cluster boundaries, may be calculated 
using ideas of two-dimensional field theory. For cluster boundaries, we show that almost 
all loops have area ClnL + O(l), where L is the size of the system, and C is a calculable 
constant. We also compute the universal ratios (A)z/(R 2 )i of the area to the squared 
radius of gyration of loops of a given large perimeter I. 

Presented at Les Houches Summer School on Fluctuating Geometries in Statistical Me- 
chanics and Field Theory, 2 Aug - 9 Sep, 1994. 



9/94 



cardy@thphys.ox.ac.uk 



Introduction. 



If we look at a typical configuration in a Monte Carlo simulation of the 2D Ising model 
at its critical point, we see a scale-invariant distribution of structures. If we focus on the 
boundaries of the clusters, there are some small ones, but others seem to stretch across the 
system. What is the statistics of the size and shape of these objects? This is a geometric 
question, not simply related to the conventional correlation functions of the Ising model. 
In fact, if we study typical configurations of the same model at temperatures T > T c , 
although it exhibits a finite correlation length in the conventional sense, the geometrical 
properties of the cluster boundaries appear to be very little different. As we shall show in 
detail, this is indeed true all the way up to infinite temperature. 

We shall discuss two useful measures of the size and shape of these objects, which, 
it turns out, may be calculated exactly for large clusters. The first is the area contained 
within a cluster boundary, which in general forms a closed loop. The second is its radius 
of gyration, defined by R 2 = (l/2£ 2 )^ ri r (n — r 2 ) 2 , where the sum is over pairs of 
points on the perimeter, and I is the number of such points. As well as the Ising model 
for T = T c and T > T c , we shall also discuss the boundaries of percolation clusters at 
and near the percolation threshold p = p c , and the problem of single self-avoiding loops 
in the plane. These all turn out to be special cases of an 0(n) model, defined as follows. 
Consider a honeycomb lattice (this is for convenience in making a clean definition of cluster 
boundaries, but similar results may be obtained for the square lattice.) At each site r 
place an 0(n) spin with components s a (r), where a = (l,2,...,n), normalised so that 
Tr s a (r)sb(r') = d a bd rr '. The partition function is 



For small x we may expand this as a sum over configurations of bonds on the lattice, with 
a weight x for each bond which is present. After taking the trace, the only configurations 
which survive (see Fig. 1) consist of self-avoiding closed loops, with a factor n that arises 
on taking the trace over the last spin in each loop. Thus 
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Figure 1. Loop gas on the honeycomb lattice, and Ising clusters 
on the dual triangular lattice. 

The advantage of the honeycomb lattice is that it has coordination number three, so 
that loops never touch at the vertices. This is not of course the usual partition function for 
the lattice 0(n) nonlinear sigma model, which is Z = Tr exp(/3 ^2 s(r)-s(r')). Nevertheless, 
they agree at high temperatures, when x ~ /3, and, for the Ising case n = 1, they agree 
exactly with x = tanh (5. Unlike the non-linear model, (jl]) makes sense for all n, and it has 
a continuous transition in two dimensions for n < 2. The phase diagram is best illustrated 
by the schematic renormalisation group flows in Fig. 2. There is a conventional unstable 
fixed point at x = x c i(n), but the whole low-temperature phase is also critical in general, 
being attracted into a stable non-trivial fixed point at x = x C 2(n). For the case n = 1, 
x C 2 = 1, corresponding to a conventional T = fixed point. 

The relation of this model for n = 1 to cluster boundaries is through duality. The loops 
in ([!]) may be thought of as the domain walls of an Ising model on the dual, triangular 
lattice. When the original model is at x = x c x, the dual spins are also critical, while 
x = x C 2 = 1 corresponds to the dual spins being at infinite temperature. Since the whole 
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Figure 2. Schematic RG flows and phase diagram of the 0(n) 
model as a function of x. 

phase x > x c \ is attracted into this fixed point, so is the whole high-temperature phase 
of the dual model. Thus the universal large distance properties of its cluster boundaries 
should be the same as those at infinite T . 

This infinite temperature Ising model on the triangular lattice may also be interpreted 
as a site percolation problem, where sites are occupied (= spin up) with probability \. 
This is precisely the percolation threshold for the site percolation problem on this lattice, 
so these cluster boundaries may also be thought of as those of critical percolation clusters, 
and their universal properties should apply to critical site and bond percolation clusters 
on any lattice. 



Area of loops. 

In two dimensions, any gauge field is confining, and as a result, Wilson loops obey an 
exact area law. One way, therefore, to express the area inside a simple closed curve is to 
regard it as a Wilson loop. To be explicit, consider a unit current J M flowing around the 
loop in its direction of orientation. Then the area is 

A loop = J ^G^(r-r')J^r)J v {r')d 2 rd 2 r' (2) 

where = (A^(r)A u (r')) is the gauge field propagator. It is simplest to take the gauge 
group to be £7(1) , so that, in the gauge A x = 0, the area is 

-4ioop = ~ J J \x - x'\5(y - y')Jy(x,y)J y (x',y')d 2 rd 2 r' 

which simply corresponds to a decomposition of the interior of the loop into strips of length 
| a; — x'\ and width 5y = by' . Although this is written in a continuum notation, in fact 
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this is valid even on the lattice if we imagine each link slightly broadened so as to form a 
narrow strip, with a current density J M such that there is unit flux along each link. This 
makes the transition to the continuum limit particularly clean and unambiguous. 

What is this current J M ? This is most easily understood by thinking of the 0(n) 
model as a complex 0(n/2) model, with complex spins S a = s a + is a+ i, and partition 
function Z = Tr + xS* ■ S + c.c). J M is then just the U(l) current, proportional to 
(l/2i)(S* ■ dS — c.c), whose integral generates the symmetry S a (r) — > e t6 S a {r). If we now 
sum over all loops and take the expectation value, we get 

n(A) = — f f \x — x'\(J y (x,y)Jy(x',y'))d 2 rd 2 r' 



2 „ 

where the (J J) correlation function is evaluated in the 0(n) model, and A is the total area 
of all loops. Note that contributions to (J ll {r)J v {r 1 )) cancel in the sum over orientations 
when r and r' are on different loops. If we use translational invariance, we finally arrive at 

(A) = ~a[ \x\(J y (x,0)J y (0,0))dx (3) 

where A is the total of the lattice. 

At its critical points, the 0(n) model, on distance scales much larger than the lattice 
spacing a is supposed to be described by a conformal field theory. J M is conserved current 
in this theory with scaling dimension 1, and, by scale invariance, rotational invariance, and 
conservation, its two-point function must have the form 

(JMMo)) = Hn) ( vvzirV ) (4) 

which is unique apart from the universal number k which characterises the theory (and the 
choice of U(l) current.) [Note that k is sometimes called the chiral anomaly, because if we 
define the left- and right-moving components of the current to be Jl,r = <h ± iJi, strict 
current conservation requires the existence of a contact term Jl{^)Jr{t') oc k5^ 2 \r — r').] 
In our case we see that the correlation function in (|3|) behaves like 1/x 2 so the integral is 
logarithmically divergent at both short and large distances. Short distances are presumably 
cut off at the order of the lattice spacing, since the continuum description breaks down 
there, but the infrared divergence is real, and indicates that the average area has no 
thermodynamic limit. Instead, if we evaluate this in a finite system of linear size ~ L, we 
should find 

(A)/A~ (fc(n)/2ra) ln(L/a) 
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The right hand side is, of course, much greater than unity for large L. This can happen 
because there are loops with loops, and thus some portions of the total area are counted 
many times in (A). 

In the same way, higher moments of A are given by integrals over higher order correla- 
tions of the current J. Fortunately, these are simple to write down, since it is an elementary 
fact of conformal field theory that the connected pieces of such correlation functions all 
vanish. The way to prove this is to consider the holomorphic and antiholomorphic com- 
ponents J = J\ — i J 2 and J = J\ + iJ 2 - Current conservation implies that a correlation 
function 

(J(z 1 )J(z 2 )...J(C 1 )J(C 2 )...} 

depends only holomorphically on the {zj} and antiholomorphically on the {Cj}), so that it 
is entirely determined by its singularities as these points approach each other. These are 
given by the operator product expansion 

J(zi) J(z2) ~ 7 — — - — ttt + regular terms 

{zi - z 2 y 

with similar relations for J. The important point is that, for a U(l) current, there is no 
term 0((zi — Z2) _1 ). The consequence is that, if we subtract off the disconnected part of 
the correlation function, the remainder has no singularities, and falls off sufficiently fast 
at infinity that it must vanish identically. 

Thus, apart from the value of k, the correlations of a U(l) current are the same in 
any conformal field theory. We are therefore free to use a free field representation, which, 
for convenience, we take to be J M = y/ke^d^rf 1 , where <f>(r) is a free scalar field with 
action Sg = (l/2*7r) §(d(f)) 2 d 2 r. Since the loops do not cross the boundary, the normal 
component of J should vanish there, corresponding to choosing Dirichlet conditions = 
on the boundary. It is also more convenient to write (||) in a covariant gauge, for which 
d^G^ = and d 2 G VLV = — £( 2 )(r — r'). It is then readily checked, using Stokes' theorem, 
that (0) correctly gives the area. 

To generate the higher cumulants of A it is convenient to consider its generating 
function, given in the free field representation by 

l -uA\ _ / —u(k/n) J J G^(r-r')e tiX e vrT dxd' rT (j}(r)4>(r')d 2 rd 2 r\ 
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Integrating by parts and using the above properties of then gives a mass term for 
the 4> field, proportional to (k/n)u, so that the result is simply a ratio of two Gaussian 
integrals 

_ uA _ det(-d 2 + (2irku/n))-^ _ -*£ m in(i+fl*^) 
{G } ~ det(-a2)-V2 " e {5 > 

where the X m /L 2 are the eigenvalues of — (laplacian), with Dirichlet boundary conditions. 
If we expand ([5]) to first order in u, we recover our previous result 



(A) = (nk/n)L 2 J2l L 



m 



since the sum has a divergence like (l/27r) 2 J d 2 p/p 2 ~ (l/27r) ln(L/a). However, the 
higher cumulants behave as 



Am 

m 



and this sum is convergent (but dependent on the shape of the region) for p > 1. 

This result implies that A/.A has mean value O(lnL), but with fluctuations 0(1). 
However, we are more interested in the area per loop, given by A/Ni, where AT/ is the 
number of loops. Since this is conjugate to the fugacity n, its cumulants are given by 



(N l ) = (nd/dn) lnZ 0(n) 

(n) 



(N 2 )-(N l ) 2 = (nd/dn) 2 lnZ 



and so on. Each of these is an extensive quantity, proportional to A, even at the critical 
point. (In fact, they are calculable exactly thanks to a series of mappings to the 6-vertex 
model.) Hence we see that Ni is of order L 2 ~ A, with fluctuations of order L. Therefore 
in the ratio A/Ni the numerator fluctuates much more strongly than the denominator, and 
it is permissible to replace (A/Ni) by (A)/{Ni). 

The final result is that, not only is average size of a loop large, but, with probability 
one as L —>■ oo, all loops have size C In L+0(1). The amplitude C is a ratio of the universal 
number k(n)/2n and the non-universal, lattice-dependent, quantity (N{)/A. From (|5p one 
may also recover the whole probability distribution P(A) of A, in the limit when o and L 
are large, by inverse Laplace transform. In particular, for 

A - (A) > L 2 , 
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we find an exponential decay 

P(A) ~ e -[*o/(27rk/n)](A/L 2 ) 

where Xq/L 2 is the lowest eigenvalue of — d 2 . 

The above results for the average area may be generalised away from the point where 
the loops are critical. In the percolation problem, this would correspond to p ^ p c , and, 
for Ising clusters, to adding an external magnetic field H. In either case the integral in (^) 
will now behave like (k/2n) ln£, where the correlation length £ behaves as £ ~ (p c — p)~^ 
for percolation, and for the Ising model as (T = T c ) and (T > T c ). 



Finite-length loops. 

Although the above calculation told us about the average statistics of loop areas, it 
did not tell us about the size or shape of individual loops. If we pick out one loop and 
measure its length to be £, how does its area and radius of gyration depend on £1 

To study this question we need a generalisation of the 0(n) model, in which the indices 
a run from 1 to n + n'\ 

(n n-\-n' \ 

1 + x^2s a (r)s a (r') + x' ^ s a (r)s a (r') 
1 n+l J 

This gives us an expansion in powers of x and x' which is a sum over configurations 
with two different kinds of mutually non-overlapping loops. As n' — > 0, the term 0(n') 
corresponds to those configurations with single primed loop: 




which defines the coefficient pg. (As n —>■ this is just the total number of self-avoiding 
loops of length £, weighted equally.) This model has a critical point at x = x' = x c . 
Turning on the perturbation (x' — x c ) breaks the symmetry down to 0(n) x 0(n'), and 
we expect the degrees of freedom with n < a < n' to become massive. 
To measure the area of the chosen loop we may use the current 

n+n' 

j' ~ (i/jm) « ds - - c - c -) » 

n+l 
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whose 2-point function is given by {J' J') = (n' /n)(JJ). Then, as before, on taking the 
expectation value in this ensemble, 

n'Y,PM)i^ = -\ [ \x\(J' y (x,0)J^0,0))dx (6) 

where {A)i is the mean area of loops of length I when embedded in the gas of critical loops 
of the 0{n) model. The integrand once again behaves as 1/x 2 , but it is now cut off at 
large distances by the length scale £' which is essentially the inverse mass of the massive 
modes. This will diverge £' ~ (x c — x') v , where v' is related to the scaling 

dimension X' of the perturbation by v' = 2/(2 — X'). Thus the right hand side of @ 
behaves like ln(x c — x') with a calculable coefficient. From this we discover the large £ 
behaviour of the coefficient of x on the left hand side: 

pi(A)t ~ (A;(n)/2n)z/ / (n)r 1 x c - 1 

Now YliPe x>£ itself is like a free energy, which, by hyperscaling, should have a singular 
part which behaves like £' 2 ~ (x c — x') 2v , which gives the large i behaviour of pg. Putting 
these results together, we find that pe ~ £-l-2v x -l and ( A } £ „ _ 

Radius of gyration. 

The radius of gyration of a single loop is given by R 2 = (l/2£ 2 ) J2 ri r 2 ( ri — r 2) 2 - We 
need to find a way of counting all loops going through two chosen bonds at r\ and r^- Let 
us define the energy density of the perturbation of the 0(n + n') model on the bond at 
r as E'(r) = Yln+i Sa ' Sa - Then E'(ri)E'(r2) has only contributions from the selected 
loops. In the limit n' — > 0, r\ and r<i must be on the same loop, since there is only one. 
Hence Yl ri r 2 ( ri ~ r 2) 2 E' {ri)E' {1-2) , when inserted into a correlation function, gives the 
radius of gyration (squared) of the loop. On averaging, we thus find that 

2n'A^2pd 2 (R 2 )ix' l ~ 2 = ^ ((n - r 1 ) 2 E'(r 1 )E'(r 2 )) ~ A ! r 2 (E'(r)E'(0))d 2 r 

where this last correlation function is to be evaluated in the perturbed 0(n + n') model. 
This last integral may be evaluated by a corollary of Zamolodchikov's c-theorem. [|1] 

The addition of the term ^2 r (x' — x c )E'(r) to the action takes the theory away from 
the confer mal point, and gives rise to a non-zero trace = of the stress tensor. In fact, 
G = 2-kv' 1 (x / — x c )E'(r), which may be demonstrated by considering the response of the 
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theory to a dilatation. |l| Thus we need the second moment of 2-point function (00). A 
simple way of deriving this is as follows: in two dimensions, conservation and rotational 
symmetry force two-point function of the stress tensor to have the form 

OV(r)T A(J (r')> = (0„0„ - 9^d 2 )(d' x d' a - gx,d' 2 )G(r - r') 

where G is a scalar function. At short distances, the theory is equivalent to some conformal 
field theory with central charge cuv-i so that 

(T„(z)T„W) ~ ~ &ff 2 G 

so G(r) ~ (cf/y/6)lnr as r — ► 0. Similarly, as r — ► oo, G(r) ~ (cm/6)hir. Now 
J r 2 (O(r)O(0))ei 2 r = J r d G{r)d r. Integrating by parts twice kills the factor r 2 , and 
the result is the integral of a total derivative, which is not zero because of the surface terms 
at r = and infinity. The final result is that 

cuv -cir = ^J r 2 (Q(r)e(0))d 2 r (7) 

In our case, the UV theory is the 0(n + n') model. In the IR, we are left with the 
massless modes of the 0(n) model. Hence the left hand side of (0) is c(n + n') — c(n) ~ 
n'(dc/dn), so that 

2j2PfZ 2 (R 2 }e x,£ ~ 2 ~ (2ttz/"V -x c ))~ 2 (4n/3)(dc/dn) 

t 

This tells us about the large i behaviour of pi{R 2 )i, and, in particular that (R 2 )e ~ ^ 2v ■ 
This implies that a loop of finite (but large) length, immersed in the gas of other critical 
loops, has a fractal dimension of v' . The amplitudes in the mean area and R 2 are not 
universal, but the ratio, in which all dimensionful and cut-off dependent quantities cancel, 
is, and we find, eliminating pi, 

(A)e k(n)/n 
(R 2 ) t ~ (dc/dn)u'(n) 

The remarkable thing about this equation is that it relates expectation values at finite 
£, that is, in the non-critical 0(n) x 0(n') model, to quantities defined at the conformal 
points. 



9 



Calculation of k, c and v' . 

It turns out that these universal quantities may be computed using Coulomb gas 
techniques. These are explained at length in the review article by Nienhuis|§ (see also 
0]), so they will be only summarised here. We first rewrite the partition function ([!]) of 
the 0(n) model in a form where the Boltzmann weights are local. This may be done by 
assigning to each vertex of each oriented loop a factor e lx or e~ lx according to whether 
the oriented walk turns to the left or the right as it passes through the vertex. In this way, 
each anticlockwise loop will accumulate a factor of e 6lx , and each clockwise loop e~ 6lx . 
On summing over orientations, we may reproduce the factor n for each loop if we choose 
X such that n = 2 cos 6x- 

The next step is to write the loop gas as a solid-on-solid (SOS) model. Define a height 
variable 4>(r), conventionally normalised so that (p(r)/n is an integer, at each site of the 
dual triangular lattice. There is 1-1 mapping between configurations of oriented loops and 
SOS heights as follows: assign <fi = on the boundary, and increase (decrease) cf> by tt each 
time a loop is crossed which goes to the left (right). The fact that all loops are closed 
makes this a consistent procedure. The resulting Boltzmann weight for each configuration 
is now a product of weights for each elementary triangle, each of which depends only on 
the differences of the 4>s at the vertices of the triangle. In the large distance limit, we now 
assume that this model renormalises onto one where the discreteness of the <ps is irrelevant. 
This may be justified a posteriori. Since the continuum limit must be a conformal field 
theory, depending on a single scalar field <fi, the simplest guess is the Gaussian model 
with action S = (g/Aiz) J (d(p) 2 d 2 r. (Note that this is different from that introduced to 
give the free field representation of J p .) The only difficult part of this argument is in 
assigning the value of g, which determines the various scaling dimensions. The argument 
is quite involved|| and we quote only the result g = 1 — (6%/7r), which corresponds to 
n = —2 cos ivg. The correct branches are given by 1 < g < 2, corresponding to 
and < g < 1, corresponding to x = x C 2- 

The free field theory is often called a Coulomb gas, because operators e iq ^ behave 
like charges of strength q: (e 1 ^^ e~ iq ^^) ~ e ~^ 2 / 3 ^ lnr . However, the SOS model is not 
completely equivalent to a free field theory, since if we calculate (e -121 *^/ 71 ") in the SOS 
model we find that it is identically equal to 1. This may be seen to occur to all orders in 
the expansion in powers of x, the first few terms of which are 

, e -i2i X ct>/^\ = 1 + e e +e e + . . . 
v 1 1 + e 6 * + e" 6 ^ + . . . 
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Since the only non-zero expectation values in the free field theory are those with total 
charge zero (as a consequence of the U(l) symmetry <p — > <p + const.), there must be a 
charge +12x/tt distributed on the boundary. The only non-zero expectation values are 
those of products of operators with total charge — 12x/tt in the interior. 

How do we identify the current of the 0{n) model in the Coulomb gas language? 
The simplest guess is to take J M = (l/7r)e Ml/ A I/ ^>, where Au is a lattice derivative, since 
this will ensure a unit current to the left or right according to whether <p steps up or down 
by 7r. However, this has the wrong charge, and in fact is not correct, as may be seen by 
computing its expectation value for a given loop, on summing over orientations: 

(J) = 1 • e 6lx + (-1) • e~ 6lx ^ 

Instead consider J M oc e At!y A I/ (e _12 * x ^/ 7r ), which does have the correct total charge. Now 

(J) oc (e~ 12lx - l)e 6lx + (e 12lx - l)e~ 6lx = (8) 

as required. However, the 2-point function cannot be (J J), since, once again, this has the 
wrong total charge. In fact, what works is (J(r)J(r')). When r and r' are on different loops 
this vanishes, because we get the same factor as in (H) on summing over the orientations of 
the loop passing through r. When they are on the same loop, we get something proportional 
to 

1 • (e~ 12ix - 1) • e 6ix + (-1) ■ (e 12ix - 1) • e~ Ux = -4zsin6x 

which may be absorbed into the normalisation of J. This is an example of how the mapping 
to the Coulomb gas is not at the level of operators, but of correlation functions. It is now 
a simple matter to evaluate (J(r)J(r')) in the Gaussian model. It is of the required form 
(§, with 

kin) = 1 ~ 9 
n 2n 2 gsm7Tg 

The central charge c(n) of the 0(n) model is known by a variety of methods. A simple 
way[|l] is to note that it is related to the free energy of a system when placed in a box of 
linear size L: 

F = — In Z ~ — (cxe /6) In L 

where xe is the Euler character. For the free field theory c = 1, but the charge on 
the boundary gives an additional contribution to the electrostatic energy. Choosing for 
convenience a disc with xe = 1, we then find 

c 1 (6 X ) 2 , 6(g -l) 2 

- = or c = 1 

6 6 g g 
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The final ingredient is z/, related to the scaling dimension X' of the perturbation 
Xm+i s a( r )s a ( r ')- In the context of spin systems, this is a symmetry-breaking pertur- 
bation inducing a crossover to lower spin dimensionality, and its analysis is somewhat 
subtle. It is first necessary to decompose the perturbation into irreducible representations 
of 0(n + n'), by writing it as 

n+n' / } \ ' 

s a (r)s a (r') - -—5 aa E(r) + -^E(r) 
z — ' V n + n J n + n 

71+1 X ' 

where E(r) = Yli +n s a(r)s a (r') is the energy density of the 0(n + n') model. The second 
term contributes merely to a shift in the critical point x c when we turn on the perturbation, 
and is not important. In the first term it is permissible to set r' = r (since the terms have 
the same symmetry) and we notice that each term in the sum transforms according to the 
representation of traceless rank 2 tensors. The scaling dimension is therefore independent 
of n', and we may use Nienhuis' result 0] for the uniaxial case n' = 1, obtained by Coulomb 
gas methods. The result is v' = 2g/(l + 2g). Finally, putting all the pieces together, we 
find 

(A) i l + 2g 
(H»)t ~ 2(1+0)* 

This result tells us about the compactness of the loops. If all loops were circular, the ratio 
would be 7r. Compact loops would have a ratio larger than this. In fact, we see for all g 
that this ratio lies between |7r and tt. For example, for Ising clusters avove T c (g = |) 
we find y^7t, for Ising clusters at T c (g = |) the slightly greater ratio ^jirr, and for single 
self-avoiding loops 0] (g = |), the value ^7r ~ 2.513. This latter value has been checked 
numerically by extrapolating the results of exact enumerations for finite I. The result is 
2.515 ... [Q, which agrees well with the analytic result. 
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